Depinning with Dynamic Stress Overshoots: Mean Field Theory 
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An infinite-range model of an elastic manifold pulled through a random potential by an applied 
force F is analyzed focusing on inertial effects. When the inertial parameter, M, is small, there is a 
continuous depinning transition from a small-F static phase to a large-F moving phase. When M 
is increased to M c , a novel tricritical point occurs. For M > M c , the depinning transition becomes 
discontinuous with hysteresis. Yet, the distribution of discrete "avalanche" -like events as the force 
is increased in the static phase for M > M c has an unusual mixture of first-order-like and critical 
features. 

PACS Numbers: 68.35. Rh, 83.50.Tq 
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A wide variety of driven systems are characterized 
by some type of elastic, extended object that is pulled 
through a quenched random medium by a uniform ap- 
plied force. If the dynamics is dissipative- e.g. sliding 
charge density waves [1], vortex lattices in superconduc- 
tors [2] and domain walls in ferromagnets [3] -such sys- 
tems exhibit a critical depinning transition in the ab- 
sence of thermal fluctuations. A small applied force, F, 
is not enough to overcome the random pinning forces 
and the system remains trapped in one of many possible 
metastable configurations. However, as F is slowly in- 
creased, some sections become unstable and move, only 
to be stopped by the elastic forces from more strongly 
pinned neighboring regions. As F is increased even fur- 
ther, there will be a sequence of these discrete, localized 
"avalanche" events with a distribution of sizes, s [4]. As 
a critical force F c is approached from below, arbitrarily 
large avalanches can occur. Above F c , the driving force 
is able to overcome the pinning, giving rise to a non- 
zero average velocity, v. The transition between the two 
phases is continuous with v playing the role of an order 
parameter that vanishes as v ~ (F — F C )P as F c is ap- 
proached from above. This class of dissipative systems 
has been analyzed by renormalization group (RG) meth- 
ods [1] with the results supported by numerical [5] and 
experimental [6] evidence, although quantitative experi- 
mental tests have been very limited. 

In some systems, in particular crack fronts in brittle 
materials [7] and motion of contact lines of droplets on 
dirty or rough solid surfaces [8], the dynamics are not 
dissipative and inertial effects can be important. In this 
Letter, we explore the effects of inertial stress transfer in 
which the motion of one segment creates a transient stress 
(in addition to the static elastic stress) on other segments. 
In the presence of pinning, the motion will be mostly 
forward; this makes the positive (forward pulling) parts 
of the transient stress transfer the most important; we 
call these stress overshoots. Although we will focus only 
on such dynamic stress transfer effects here, the more 
direct effects of inertia will have similar consequences: if 
a segment moves forward in an undcrdamped manner so 



that it overshoots - passing a new local minimum before 
relaxing back into it - the stress it transfers to the other 
regions will reflect this overshoot. 

Our focus will be on the nature of the depinning tran- 
sition in the presence of stress overshoots; a key ques- 
tion is whether or not the transition from the pinned to 
the moving phase remains continuous or becomes "first- 
order"- characterized by a finite jump in v as a function 
of the applied force and perhaps by hysteresis. We will 
also consider the distribution of avalanche sizes as the 
depinning transition is approached from below. 

Near the depinning transition, the motion will be very 
jerky with a lot of starting and stopping. Thus, the es- 
sential physics can be captured by models in which space, 
time and the manifold's displacement in the direction of 
motion, /i(x, t), are all discrete. The elastic stress on a 
segment x takes the general form 

<r(x, t) = J2 ^( T ) ft (y> t-r)-j h(x, t), (i) 

y t> o 

with the static stress transfer given by J xy = ^ T Jxy(r) 
and J — J2 y ■ The simplest form of the stress over- 
shoots is when they only last for one time step, i.e. 
Jxy{T > 2) = 0, and J xy is non-zero only for nearest 
neighbors. In this case we write, </ xy (0) = ^ 1+ z^ and 
J xy (l) = — 4|- for x/y where Z is the number of near- 
est neighbors and M is the magnitude of the stress over- 
shoot. With this dynamic stress transfer, the motion of 
any neighboring segment, y, causes an extra stress on the 
xt/i segment for one time step only. For M < 0, a key 
property of the elastic stress transfer is that it is mono- 
tonic, i.e. increasing in time when the /i's increase. This 
monotonicity property can be used to prove the unique- 
ness of the average velocity and thereby a unique critical 
force [9]. It is also an essential feature for the validity 
of the RG analysis [1]. For M > 0, the stress transfer 
is non-monotonic in time so that the average velocity is 
not necessarily a unique function of F. 

We take the displacements /i(x, t) to be pinned at a dis- 
crete set of possible values with associated pinning (yield) 
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strengths fy [x, ft.(x, t)]. The dynamics is simple: if the 
total force on segment x exceeds fy at that point, the 
segment jumps forward, i.e. if cr(x, t) + F > fy[x, /i(x, t)] 
then h(x,t + 1) = h(x,t) + A[x, h(x,t)} with A the dis- 
tance to the next pinning position; otherwise /i(x, t+l) = 
h(x,t). As long as the jumps A[x, h(x, t)} are random 
variables independently drawn from a smooth distribu- 
tion D(A)<iA, the randomness in the pinning strengths is 
not essential and, for simplicity, we make them all equal. 

We now take the first step in analyzing such systems by 
studying a mean field- more precisely an infinite-range 
-model. For the case without stress overshoots, such a 
model was the needed starting point for the RG analysis 
of finite-dimensional systems. We thus consider a model 
of N segments each coupled to all the others. The stress 
on h(x,t) is then simply 



a(x, t) = h(t) - h(x, t) + M (h(t) - h(t - 1)) 



(2) 



where h(t) = ^ J2 y h(y,t) is the average over the other 
segments. 

In this infinite range model, all sites are statistically 
equivalent and coupled only via the mean field, h(t). 
We can therefore fully characterize the configuration 
by the distribution of the excess force on a segment 
ft(x) = a(x,t) + F — fy- Dropping the x index, we 
denote this distribution, p t (ft)df t . Since each segment is 
either pinned, ft < 0, or will jump forward by some ran- 
dom amount, A, the equation of motion for the excess 
force distribution is 

/>oo 

Pt+i(ft+i)= / d(/> t 5(<pt+vt+i+M(vt+i-vt)-ft+i) 



x{pt(&)e(-&)+ / D(A) Pt (cp t + A)e^ t + A)dA}, (3) 
Jo 

where <pt is the new excess force on a segment before the 
others have moved, vt = hit) — hit — 1) and 8 is the unit 
step function. The self-consistency relation for vt follows 
from the "spatial" averaging of the excess force: 



<ft >-- 



ftPt(ft)df t =F - fy + Mvt. (4) 



For static solutions, v = 0, the stress overshoot plays no 
role and one can show that there are many solutions as 
long as the applied force is less than the static critical 
force F c q — fy — where the bars here denote aver- 

ages over D(A). We will return later to the response to 
an increase in the force in this regime. For now, we turn 
to the moving phase. 

We first consider the steady state limit for which 
v t = v > is independent of time. Equation (^) im- 
plies that in this case the stress overshoot has the same 
effect as an additional applied force, F — > F + Mv, and 
the steady state p(f) depends only on v (and D). It 
is easiest to work with the particular jump distribution 



D(A) = e~ A which corresponds to random pinning posi- 
tions with density 1/A chosen to be unity. However, the 
qualitative properties of the v(F) curve are applicable 
to any well-behaved, smooth pinning distribution with 
-D(O) finite and non-zero. For M < 1 and F > F c0 , v 
rises continuously as 



F — F, 



cO 



1 - M 



(5) 



for v << A = 1. The fact that v is linearly proportional 
to F — F c in mean field models with jumps has been 
derived previously for dissipation-dominated systems [1]; 
it is due to the fact that above the depinning transition 
the average displacement of a segment responds linearly 
in a non-singular way to an increase in the total force on 
it. Note that as the applied force is increased further, the 
v of the model is constrained in an unphysical way by the 
average distance between pinning sites. As a result, v will 
saturate in a manner that depends on the particulars of 
the pinning distribution. For the exponentially decaying 
jump distribution, the full expression for v is 



-e + (i + ^(ep + /i) 2 + 2e 
1 + 2/x 



(6) 



with e = F — F c0 and // = M — M c . 

As M increases, the amount of applied force required 
to reach a specific v decreases. When M reaches a critical 
value of M c — 1, the width of the continuous depinning 
transition shrinks to zero. When both the parameters M 
and F are critical there is a "tricritical" point. Near this 
tricritical point the scaling behavior is v ~ \fi\ ~ y^O' m " 
dicating a new universality class of depinning transitions. 

For M > M c , there is a jump in v from v m i n ~ /x to zero 
when F is decreased through a new lower critical force 
Fj — F c0 — (1 + /j, — yjl + 2(i). The two stable solutions 
in this regime consist of v = and a branch with jp < 0. 
There is also a third solution, but it is unstable. Between 
Fj and F c q there is hysteresis as might be expected at 
a first-order like phase transition. A greater decrease in 
the force is required to stop the motion in the presence 
of overshoots than in their absence, but once the overall 
motion does stop the overshoots will have less effect as 
F is increased back up again. 

We now investigate the avalanche dynamics below the 
depinning threshold. An avalanche is the forward motion 
of a finite number of segments in response to the initial 
motion of one segment. The applied force is increased 
only in order to trigger the motion of the initial segment 
and is then held fixed until the avalanche stops. The 
segments that move forward in a single finite avalanche 
are those whose / is within an infinitesimal- O(i)- slice 
near / = of the distribution p{f); for continuous 
all that matters is thus p(Q). The total number of seg- 
ments that hop forward at one time step, n*, will have 
poissonian statistics whose mean is determined by the in- 
crease in the total force on a segment from the previous 
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time step. The mean number that hop at a given time t 
is given by 

<n t >=p(0)A[(l|M)n t _i-Mn f _ 2 ]. (7) 

If this mean becomes negative, the avalanche will stop. 
We define the parameter, T = p(0)A, which is essentially 
the mean local response to an increase in the total force. 

We are primarily interested in the large avalanches, 
for which n t will typically be large. Given large n t -2 and 
n t _i, the randomness will cause approximately gaussian 
fluctuations in n t of magnitude of order y/n~t- With no 
stress overshoots, M = 0, we can approximate the dy- 
namics in continuous time in the Ito representation as 



-A2~(r-l)n(t) + v ^(t)»7(t) (8) 

with r] white noise. The constraint that n(t) be posi- 
tive makes finding the distribution of the avalanche size, 
s = J dtnt, a first passage problem. A critical point oc- 
curs at T = 1 (p c (0) = ^), at which there is a power 
law distribution of avalanche sizes, Prob(ds) ~ s-ds for 

large s [10]. Beyond this point, infinite avalanches occur 
as the motion of one segment spawns more than one de- 
scendant at the next time step. For smaller T, a cutoff 
in the avalanche size distribution given by s <~ (l-r)' 2 ( 
which is not equal to the average as the distribution is 
very broad), and Prob(ds) <~ -^-e - " r~ ds for large s 
as shown in Fig. 1. 

In the more general case with stress overshoots, the 
stochastic equation becomes second-order in time since 
the number of segments that hop forward is dependent 
on two previous time steps and we have, approximately, 

rM^, ( r-i W t) + (rM-i)M 

+ ^W)v{t) • (9) 

When we introduce a small amount of stress overshoot, 
the avalanche dynamics is strongly overdamped. The 
critical point is still at T = 1 but the distribution of 
avalanche sizes is modified by the coefficient (M — 1). 

When the stress overshoot magnitude is increased to 
the critical value, M c = 1, the avalanche dynamics be- 
comes undamped resulting in a new exponent for the dis- 
tribution of avalanche sizes. At this tricritical point, wc 
have d — y // n(t)r](t) which maps the distribution of 
avalanche sizes into a first passage problem of a particle 
with random acceleration. The distribution of avalanche 
sizes broadens to Prob(ds) <~ -\ds for large s [11]. The 
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typical avalanche of size s now lasts for a time r ~ in 
contrast to t ~ s5 in the absence of stress overshoots. 
The scaling form of the distribution within the vicinity 
of this tricritical point is 



Prob(ds) ~ -^j l F{s{M - M c ), s 1/2 (l - T))ds. (10) 

For M > M c , the deterministic part of the avalanche 
equation of motion is that of a time- reversed harmonic os- 
cillator. If p(Q) is small- i.e., below the critically damped 
line —n(t) grows exponentially with time before crash- 
ing down to zero as a result of the oscillation. Above 
some critical T, T C (M) = j^jji (from Eq. (7)), large 
avalanches grow exponentially and become infinite. On 
the critical line, some avalanches become infinite while 
others remain finite. Because of the exponential growth 
of n(t) both above and below the critical line, the ran- 
domness will not have much of an effect once an avalanche 
becomes large and a deterministic analysis where all the 
randomness is encoded into random initial n and ^ 

at 

becomes valid. Just below the critical line, there is a 
clear delineation in the distribution of avalanche sizes 
between those avalanches that would remain finite and 
those that would become infinite at the critical line. In- 
deed, just below r c (M), a finite fraction of the avalanches 
will have size of order § ~ exp(A(M)/ y/(T c (M) - F)) to 
within a multiplicative factor of order unity; these last 
roughly half the oscillation period. Note that for M > 1, 
A ~ ji = M — 1 and the scaling of 1 — T versus p is like 
that of F — F c0 versus p in the moving phase. Most of 
the other avalanches will be small, but with a long tail 
in log s scaling as 

P TO b(ds)~-^ds=f^l (11) 
s(ms)^ (ms) z 

for 1 << Ins << Ins. This is an unusual distribution 
with In s being the natural variable. This broad distri- 
bution for the finite avalanches persists at the depinning 
transition, but now there is also a non-zero probability 
that an avalanche will be infinite as indicated in Fig. 
1. These results suggest a new hybrid transition for 
M > M c with a diverging scale, s, but also disconti- 
nuities, particularly in v(F). 

So far, we have not established the connection between 
the driving force, F, and the density of about-to-jump 
segments p(0). This is needed to connect together the 
analysis of the pinned and the moving phases. As long as 
the avalanches remain finite and the driving force is adi- 
abatically increased, the total force distribution p(f)df 
will be independent of M; the avalanche distribution will 
only affect when segments jump and hence the small scale 
{f ~ jj) behavior of the stationary configurations be- 
tween avalanches. For generic initial conditions — con- 
tinuous p(f) which must be consistent with randomly po- 
sitioned pinning values of h(x, t) — T will increase with 
F to reach unity only right at F c0 , although the way in 
which it approaches this value is non-universal and his- 
tory dependent. Thus for M < 1, T c — 1 and the infinite 
avalanches will only occur when F > F c q, the same point 
above which v > 0. 
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FIG. 1. Log- log plot of the probability of an avalanche oc- 
curring of size within the interval [s/y/2, \/2s). The circles 
are numerical data in the absence of stress overshoots, with 
solid circles at the critical T = 1 and open circles just below, 
r = 0.999. The squares are at the critical M = M c — 1. 
The closed squares are at the tricritical point r = 1 with the 
open squares just below, T = 0.999. The solid lines have slope 
^ and j, representing theoretical predictions for the critical 
M < 1 and tricritical M = 1 points. The open diamonds are 
for M = 2 and T = 0.88, just below the "first-order" line; 
these should roughly match the dashed curve which is from 
the deterministic approximation with an exponential distri- 
bution of initial velocities. The dot-dashed curve represents 
Eq. (11). 

For M > 1, however, T will reach the value Y C (M) at a 
value of the driving force Fj , that depends on the initial 
conditions as well as on M. From an analysis of Eq(3), it 
can be shown that Fj is between F*(M), the minimum 
F for which a moving solution exists, and F c q, as should 
be expected. Thus, if F is increased above Fj and back 
down there will be hysteresis. If F is decreased below F£ 
so that the system stops and is then increased back up 
again, the behavior can be complicated as p(f) will gen- 
erally not be continuous. Depending on the parameters 
(including D(A)) and the history, the depinning transi- 
tion on the second increase of F can either be of the hy- 
brid type discussed above with an exponentially growing 
characteristic size, s, of large avalanches which diverges 
at a new F} ; or a strongly discontinuous transition with 
r jumping from below r c to above discontinuously and 
no extensive number of precusor large avalanches. 

We have found that in a simple infinite-range model of 
depinning with dynamic stress overshoots, rather sub- 
tle behavior can occur. For large overshoots, a hy- 
brid transition exists when the force is increased adia- 
batically. This is characterized by exponentially-large, 
almost-deterministic avalanches followed by runaway at 
a critical force to a non-zero velocity state; this moving 
state persists to a lower critical force as the force is de- 
creased back down. Various history dependent effects can 
occur in this large overshoot regime. For smaller stress 
overshoots, there is a reversible critical transition similar 



to that in the absence of overshoots. These two regimes 
are separated by a "tricritical" point which represents a 
new universality class. 

All of this behavior should obtain for a broad class of 
infinite-range models, including ones with locally under- 
damped relaxation caused by inertia [12,13]. But the cru- 
cial question is: what of this behavior persists in finite- 
dimensional models? For large stress overshoots the de- 
pinning is surely of a different character than in purely 
dissipative systems, although what aspects of it might 
resemble the infinite-range model is far from clear. For 
small stress overshoots, it has been argued elsewhere that 
the critical force on increasing M is generically less than 
F c q in finite-dimensional models with long-range inter- 
actions [14]. This is the relevant situation for fronts of 
planar cracks propagating through heterogeneous brit- 
tle materials and one would guess that the same should 
hold more generally. But whether this implies that the 
transition immediately changes character when any stress 
overshoots are present - and thus presumably that there 
is no tricritical point - or whether dissipative-like critical 
behavior can still exist with small overshoots, we must 
leave as an open question. 
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